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ABSTRACT 

It has been known for some time that there are two methods to calculate Bk with staggered 
fermions: one is the two spin trace formalism and the other is the one spin trace formalism. 
Until now, the two spin trace formalism has been exclusively used for weak matrix element 
calculations with staggered fermions. Here, the one spin trace formalism to calculate Bk 
with staggered fermions is explained. It is shown that the one spin trace operators require 
additional chiral partner operators in order to keep the continuum chiral behavior. The 
renormalization of the one spin trace operators is described and compared with the two spin 
trace formalism. 
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Abstract 

It has been known for some time that there are two methods to calculate Bk with 
staggered fermions: one is the two spin trace formalism and the other is the one spin 
trace formalism. Until now, the two spin trace formalism has been exclusively used for 
weak matrix element calculations with staggered fermions. Here, the one spin trace 
formalism to calculate Bk with staggered fermions is explained. It is shown that the 
one spin trace operators require additional chiral partner operators in order to keep 
the continuum chiral behavior. The renormalization of the one spin trace operators is 
described and compared with the two spin trace formalism. 



1 Introduction 

From the Standard Model, one can derive the low energy effective Hamiltonian of electro- 
weak interactions by decoupling heavy particles such as the W^, Z bosons and the t, b, 
c quarks [|T], 0, 0, ^, |^. The resulting effective Hamiltonian is composed of four-fermion 
operators and bilinear operators. 

In order to deduce the CP violation phase in the CKM matrix from the measured weak 
parameter e, we require the Kaon B parameter [0, ||, |]. To extract Bk, we need to know 
the weak matrix element M.k defined as 

Mk = I [si, {I + l,)d][si,{^ + 75)^^] I K') ■ (1) 
Bk is the ratio of J^k over its vacuum saturation value, 

R ^ Me (2) 

l{K^ I [si.i.d] I 0)(0 I [si,i,d] I K^)- ^ ^ 

In order to calculate A4k on the lattice, we need to find a set of operators which can describe 
the same physics on the lattice as the continuum four-fermion operator. 
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Especially for the weak matrix elements involving the pseudo-Goldstone bosons (for ex- 
ample pions, kaons, etc), it has been preferred to take maximal advantage of the Ua{^) sym- 
metry of the staggered-fermion action which is not manifest in the lattice Wilson fermion 
action]^, |^ . This is the reason why we choose staggered fermions for our weak matrix element 
discussion in this article. 

There are two methods to transcribe a continuum weak matrix element (for example Bk) 
onto the lattice one is the one spin trace formalism and the other is the two spin trace 
formalism. The four fermion operators can be expressed as products of operators bilinear in 
the fermion fields. In the one spin trace formalism each external hadron is contracted with 
both bilinears of the four-fermion operators simultaneously. In the two spin trace formalism 
each external hadron is contracted with only one of the bilinears in the four-fermion operators 
P, 1^. By Fierz rearrangement, these two choices describe the same continuum physics. 

Until now, the two spin trace formalism (2TR) has been used exclusively for calculations 
of weak matrix elements with staggered fermions in lattice QCD simulations |^ 0) §] • Iii this 
article, we would like to explain the one spin trace formalism (ITR) which can also be used 
to calculate weak matrix elements (for example Bk) with staggered fermions 0, |^. 

Let us summarize the main contents in this article. There is a difficulty in transcribing 
the one spin trace operator to the lattice in a way which preserves the same chiral behavior 
and the same leading logarithmic behavior as the continuum AS* = 2 operator. The key 
point is that we must add specific operators named chiral partner operators to each channel 
in order to keep the correct chiral behavior. The question then is whether the additional 
chiral partner operators can still guarantee the same leading logarithmic behavior as the 
continuum operator. It is shown that the answer is yes. For a specific representation of the 
exact staggered fermion chiral symmetry group t/^(l), it is shown that both correct chiral 
behavior and correct leading logarithmic behavior are guaranteed. The chiral behavior and 
the one-loop renormalization of the operators in both formalisms will be compared with 
each other. In order to make the arguments more clear and specific, we will restrict our 
discussion to B^- Most of the arguments can be extended simply to other weak matrix 
element calculations in the one spin trace formalism on the lattice. 

This paper is organized as follows: in section 2, we introduce the notation for staggered 
fermions, definitions for bilinear and four-fermion operators and the Fierz transformation of 
operators on the lattice. In section 3, we explain how to transcribe operators for Bk both in 
the two spin trace formalism and in the one spin trace formalism. In section 4, we derive the 
chiral Ward identities and the chiral limit of the continuum and lattice operators. We also 
introduce the concept of chiral partner operators for the lattice operators in the one spin trace 
formalism. In section 5, we carry out the renormalization for the bilinear and four-fermion 
operators. Then we connect the lattice operators with the continuum operators at the one 
loop level. We also explain the effect of the chiral partner operators on the renormalization 
of the one spin trace operators. In the end, we give a summary and conclusions. 

2 Notations and Terminology 

In this section we will specify our notation for the action, fermion fields and composite 
operators on the lattice. Also we will introduce the concepts of one spin trace operators and 
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two spin trace operators. We will explain the difference in the Fierz transformations of the 
lattice and continuum four fermion operators. 



2.1 Action and Operators 

The action for the staggered fermions in the notation used in Ref. 
1 



IS 



2a 



[x{n)Uf,{n)x{n + fi) - x{n + fi)Ul{n)x{n)] + mx{n)x{n) 



(3) 



where n = (rii, 712, tt-s, 714) is the lattice coordinate and i]^{n) = (— l)"!"* The single- 

component staggered fermion field x corresponds to 4 degenerate fermions (c/. fermion 
doubling problem). Golterman and Smit proposed one method to express the staggered 
fermion field as 4 degenerate flavors in terms of Euclidean Dirac spin and flavor matrices 
obtained by chopping momentum space into 16 divisions [0. Kluberg-Stern, Morel, Napoly 
and Petersson suggested that we can take hypercubes in the coordinate space as the basic 
units which correspond to the individual points of a lattice two times coarser and interpret 
the 16 hypercube coordinates as 4 Dirac spin and 4 flavor components []TT|. In Ref. [ |11| , 
Kluberg-Stern et al. defined the Dirac field Q{y) as 



1 



NfJNf 



Ei^ArxivA) 



(4) 



where a is the Dirac spin index, i is the flavor index and Nf is the number of degenerate 
flavors {Nf = 4). Define also 



and 



yA = 2y + A, with A G {0, 1}^ 



lA = itltlt'lt' 



(5) 



(6) 



If we set the gauge links equal to unity, the staggered fermion action can be expressed in 
terms of the Q field as follows: 



^ = «^im[x(2/A)(7M®/)AB 

A,B 



'x{{y + 1^)3) -x{{y- i^)BY 

Aa 



E E [xiVA) (/ ® I)ABX{yB)] 

x{{y + + x{{y - P')b) - 2x{yB) ^ 



A,B 



A,B 



'Q{y + i^) -Q{y- f^y 



4a 



(2a)^ r 
] + {2aymY.[Q{y){I ^ I)Q{y)] 



(7) 



-f(2a)^E[^3(l/)(75®U) 



'Q{y + fi)+Q{y-fi)-2Q{y)\ 



i2ay 
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where fl represents a unit displacement on the coarse lattice and the matrices 'Js^^f represent 
the standard tensors (or direct products). The matrices 7s ® C,f are given by 

(75 ® ^f)ab = ^Tr(7]i757B7],) . (9) 
These two matrices are related to each other by 

{is ® ^fT''^^ = ^ E iTii^^hs' , (10) 

^ A,B 

where a, (3 are Dirac spin indices and i, j are flavor indices. This interpretation is called 
coordinate- space method. 

It is known that the flavor interpretation in the coordinate-space method by Kluberg- 
Stern et al. can be related to the momentum-space method by Golterman and Smit [12] as 
follows: 

1 , 1 



{is ® iF)AB = E Tri-^)^-''ilS ® ^Fho^i-ir^" . (11) 



The momentum-space matrices (75 ^f)ab introduced by Golterman and Smit are related 
unitarily to the coordinate-space matrices (75 <S) ^p)"^'^''^^- Here we adopt the coordinate- 
space method. 

The continuum limit of the staggered fermion action on the lattice is a kind of QCD with 
four degenerate flavors {Nf = 4) |T^. In the continuum, operators of the form 

q{x){ls ® ^F)q{x) (12) 

can be used as interpolating flelds for fermion bilinear objects 0], where 75 is the Dirac 
spin matrix and C,f belongs to SU(4) flavor symmetry group. There are a lot of choices to 
transcribe the lattice operator for a given continuum operator. The conventional choice of 
bilinear operator which we will follow here, is 

OsF{y) = -^J2x{yA){is ^F)ABl^{yA,yB)x{yB), (13) 

AB 

where lA{yA,yB) is a product of gauge links that makes Osf gauge-invariant |jl3|, 0, |I4|- 
There are two simple choices for U{yA, yB)- One is that U{yA, yB) is replaced by the identity 
and the spinors xijJA) and xijJB) are evaluated in Landau gauge (so called Landau gauge 
operator). The other is that the link matrices are included and that Osf is made as symmetric 
as possible by averaging over all of the shortest paths between ua and ys [0, 0- 

i^{yA,yB) = ^^Y.u(yA,yA + Api)--- 

4. p 

U{yA + Api + Ap2 + Ap3 , vb) , (14) 
where P is an element of the permutation group (1234) and 

= {B, - A,)[i . (15) 
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There are two kinds of four-fermion operators which have different color contraction struc- 
ture. The general form of color two-trace operators is 

[x}SyA)hs ® ^F)ABX%{yB)][x%{yc){is' ^F')cDxjSyD)]^ab{yA,yB)i^^^ (16) 

and the general form of color one-trace operators is 

[X/i(yA)(7S®^F)ABX/2(yB)][X/3(yc)(75' ®^F0cDxl(yD)]Wad(yA,yD)Wc6(yc,l/B) (17) 

where /i, /2, /s and label the continuum flavor (for example, u, d ,s, ■ ■ ■) and the two 
choices for Uab{yA,yB) are the same as given for the bilinear operators. The reader should 
recognize that at this point we have introduced (following Sharpe & Kilcup |^) flavor 
in two ways: First the 4 degenerate flavors associated with each staggered fermion fleld 
x{yA), labeled by the index A, and second the possibly non-degenerate additional flavors 
corresponding to distinct staggered flelds x/j labeled by the indices /. 

We can represent the four fermion operators in Bx in terms of S (scalar), V (vector) , T 
(tensor), A (axial) and P (pseudo-scalar) as follows: 



[X% ( 7m ^ I) X%] [X% i l, ® I)X %] 

[/I X OJ^^C2;C3C4 



JciC2;c3C4 



[X/U7.5®/)X/^][X/^(7.5®/)x}: 



(18) 
(19) 
(20) 
(21) 



: m\{i,^®^^)x%][x%h,^®^^)x%] 

where /i (1 < i < 4) represent the continuum flavor (u, d, s, ■ ■ ■), Cj (1 < i < 4) represent 
the color indices in the fundamental representation of SU(3) and for example [A x PYcJc2\c^ct 
represents a four-fermion operator on the lattice which has the axial spin structure A and 
the pseudoscalar-like flavor structure P. 



2.2 Fierz Transformation 

Let us introduce the following notation : 

8 = 101, V = ^ 7^ (g) 7^, T = ^ cr^^ (g) a^j,. 



^ = H ® 7m5, ^ = 75 ® 75 



(22) 



where ® means direct product and a^y = |[7^,7;^]. The same notation also applies to the 
flavor structure ® by switching 7 to C,. 

In terms of the notation given in Eq. (^), the Fierz transformation of S, V, T, A and 
P can be expressed as follows: 

{T,,T2,Ts,T,,T,) = {S,V,T,A,P) 

(Xi)a,/3';a',l3 = Pij (T j) a,l3;a' ,l3' (23) 
1 1-1-11 

4-2 0-2 -4 

-6 -2 -6 (24) 

-4 -2 0-2 4 

1-1-11 1 
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From this relationship, we can show that the following two operators are identical: 

S^^,{1 - ^,)D^S\il - ^,)D' = 5"7^(1 - ^,)D'S\il - ^,)D\ (25) 

The continuum A^" = 2 operator is Fierz-transformed into itself. 

Let us use the above relationship (Eq. (0)) to obtain the Fierz transformation of the 
lattice operators for staggered fermions. In terms of Q"'*, the spin matrix and flavor matrix 
can be separated completely. Therefore, the spin and flavor matrices are Fierz-transformed 
separately according to the relationship (Eq. (p3|)). The Fierz transformations for the lattice 
four-fermion operators are 

for r„r„r;,r, G {^,v^,T,AP} 

[r. X T^J^:^ = _F^^ F,, [F, X T^^:^:f^;i: (26) 

where Fj, F/ represent the spin structure and F^, F^ represent the flavor structure as explained 
in Eq. (|TB]-0). 

Now let us obtain some useful relations which will important later. From the lattice Fierz 
transformations Eq. (P^D , we obtain the following relationship on the lattice: 

i [VxS + AxP + VxP + AxS]ilit:S£ 

= ^^[{V + A)x{S-T + P)]4\t;Si (27) 

which will be used later to show that the lattice operator in the one spin trace formalism is 
different from the corresponding operator in the two spin trace formalism, as well as to show 
that the Fierz transform property of the continuum will be recovered in the continuum 
limit of a = 0. This is because the operators with different flavor structure from that (^5) of 
the pseudo-Goldstone pion, {[{V + A) x S] and [{V + A) x T]) can not contribute to in 
the continuum limit of a = 0. 



3 Operator Transcription for Bk 

In this section we explain how to choose the lattice operator in order to calculate the weak 
matrix element found in Bk- We follow quite closely the original work on the two spin trace 
formalism of transcribing operators for the weak matrix elements done in Ref. [|^, |T3|. In 



this section, we explain and compare two spin trace formalism and one spin trace formalism 
in detail. 

The continuum operator related to B^ 0, |l| is 

0^=2 = r7M(l-75)ci1[^V(l- 75)^1 (28) 
= [s-^^{l-j,)d'][s\{l-^,)d'^] . (29) 

As you can see from the above equations, the matrix element in B^ is Fierz-transformed 
into itself in the continuum. Since there are four degenerate flavors on the lattice, the 
Fierz transformation is completely different on the lattice from that in the continuum. This 
difference allows two different ways of transcribing operators on the lattice. 
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The continuum AS* = 2 operator has only two flavors (s and b). Following Ref. 0, 
we may also introduce four valence quark flavors {S, S' and D, D') to represent the same 
physics: 

{K' I [^^^7^(1 + i,)d^][s\{l + 75)^1 I K') (30) 
2{K' I [S'^7m(1 + 75)^1 [^'V(l+75)/^1 I K'') 
+2(/?° I [S^^,{1 + ^,)D%S'\{1 + 75)/^'1 I K'') (31) 

+2{K' I [S^^,{l + ^,)D"][S\{l+^,)D^] I (32) 

where the primed and unprimed flavors are supposed to have the same mass. We have 
introduced additional species S' and D' to indicate a particular contraction. Similarly, the 
hadronic eigenstate is composed of the valence quarks S and D, while K'^ contains S' 
and D'. In the Eq. (0), the first bilinear has to be contracted with K and the second 
bilinear with K'. This is the so-called two spin trace contraction. Another possibility is the 
one spin trace contraction as in the Eq. ( |32D where the external kaon is contracted with 
both bilinears. 

As you can see in the above, the Fierz transformation connects Eq.(^) with Eq.(|5^). 
They are identical to each other. But the lattice version of Eg . (|3ll) and Eg . (|32|) can not 
be related by Fierz transformation. Thus, the righthand side of Eg. ( piD corresponds to two 
spin trace formulation of the operator on the lattice, while the righthand side of Eg. (p^) 
corresponds to one spin trace formulation, as will be described later. 



Mk = 
-^Mk = 

or — > Mk = 



3.1 Two Spin Trace Formalism for Bk 

The detailed explanation of two spin trace formalism for Bk is given in Ref. 0]. 

First of all we need a Kaon operator on the lattice. Since the conserved axial current has 
flavor ^5, it is better to use the composite operator of the lattice pseudo-Goldstone boson 
with the flavor matrix ^5 as a Kaon operator: 

K{n) = Xdin)e{n)xsin) 

K{t)=Y,K{n,t) (33) 

n 

K\n) = Xd'in)e{n)xs'{n) 

K\t)^Y.K\n,t) (34) 

n 

where X^n represents the summation over spatial volume needed to construct an operator for 
zero-momentum Kaon. We need to specify the four fermion operator on the lattice which 
corresponds to the continuum four fermion operator in Eg. (|3l|) . The two spin trace formalism 
means that the operators consist of two bilinears and that each external kaon is contracted 
with one of these bilinears. Hence the four fermion operator is 

O^^l = (V X P)lll + {V X P)ZS, + (A X P)lll + (A X P)Z^, (35) 
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with 



2TR — 
ab;ba 




(36) 


2™ _ 
aa;bb — 


1 d 'd' 
jyl X P]aa-bb 


(37) 


2TR — 
ab;ba 




(38) 


2TR _ 
aa;bb 




(39) 



{A X P) 
{A X P) 

where a and b represent color indices in the fundamental representation of SU{3) and the 
notations in the righthand side are as defined in Eqs. (|TH[]^). As you can see in the above, 
we choose the vector and axial channel in such a way that the flavor structure is the same 
as that of the pseudo-Goldstone boson (Kaon) and the contraction with the external Kaons 
has a non-vanishing flavor trace. The operators with different flavor structure from the 
pseudo-Goldstone boson (Kaon) are supposed to vanish when it is contracted with external 
Goldstone Kaons, because of vanishing flavor trace. However, one needs to note that the 
argument of vanishing flavor trace holds only for tree-level discussion (i.e. in the classical 
limit), because the radiative corrections in lattice QCD mix an operator with operators of 
different flavor matrices. This point will become more clear later when we discuss about the 
renormalization. 

The two spin trace form of A4k is prescribed by 

MY'^{t,t) =< K{t)0^^lK'{t ) > . (40) 
3.2 One Spin Trace Formalism for Bk 

For the one spin trace formalism, we use the same operator for the kaon as in the two spin 
trace formalism. The 4-fermion operator in one spin trace form, however, is 

o^^'^ = iy X p)Zl + {V X p)Yj, + (A X pyji + (A X p)Z^, 

chiral partner ^ ' 

with 

{V X P) 
{V X P) 



ITR _ 
ab;ba 


2Nf ^ ^^'"'''^'^ 


(42) 


ITR _ 
aa;bb 


2Nf ^ ^^"'"'^^ 


(43) 


ITR — 
ab-.ba 


2n}-^ ^ ^^"^'^^ 


(44) 


ITR. _ 
aa\bb 


2A^/^^ ^ ^^'^"'^^ ' 


(45) 



where the notations in the righthand side of the equations are defined in Eqs. (|TB|-^T]) 
and O^^^g^j partner determined explicitly later when we discuss the chiral behavior. 
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In the one spin trace case, the external kaons have to be contracted with both bihnears 
simultaneously. The one spin trace form of Aix is defined as 



the two spin trace case, for any choice of the flavor matrix F (i.e. ^p), the contraction of the 
external Goldstone Kaon fields with the {V x F) and [A x F) does not vanish. However, 
in order to take maximal advantage of Ua{^) symmetry, we will choose the axial channel as 
{A X P). In order to achieve the same leading logarithmic behavior in the one spin trace 
form on the lattice as the continuum operator, the {A x P) operator needs {V x P) as a 
vector channel. This point will become clear later when we discuss renormalization. 

In order to achieve the same chiral behavior for the one spin trace form on the lattice as 
the continuum AS = 2 operator, it is required that once we choose the fiavor of the axial 
current, we should add a vector current which corresponds to the chiral partner of the given 
axial channel (for example, {A x P) operator needs {V x S) operator as a chiral partner). 
Once we add the chiral partner to the given axial channel, the operator in the one spin trace 
formalism can have the same chiral behavior on the lattice as in the continuum. The point 
is that for the vector and axial channel to satisfy the appropriate chiral limit, the correct 
chiral partners must be added. 

4 Chiral Behavior 

With the two lattice prescriptions for Bx introduced in the previous section, we have to prove 
that each formalism separately possesses the correct chiral behavior. In the continuum, one 
can show that / dt dt A4x{t,t ) vanishes both in the limit of = rris and the limit of 
m^i = iTigi . For our proof, we will have to use several lattice Ward identities which will be 
explained in the following sections. 

4.1 Continuum Chiral Ward Identity 

In the continuum we have the following Ward identity 




(m, + m,) j d%{0 I T[S''{x)S^iy)i^,)p,^D'^iy)D^iz)] \ 0) 
= (75).p(0 I T[D^{x)DP{z)] I 0) + (0 I T[S^{x)S''{z)] \ 0){^,),^ , 



(47) 



where T indicates the time ordered product. 

Using Eg. (|47|) , we can obtain the following relationship 




(48) 



C=V,A 



where 



]Cont 
AS=2 



(49) 
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with 



^ J d'yD{y, th,S{y, t) (50) 

K'Cont^^ = / £yD>^y^ t)j,S'{y, t) . (51) 



The above result corresponds to the chiral hmit of the two spin trace formahsm. Applying 
a Fierz transformation on the right-hand side of Eq. (^8]), we obtain another form of the 
Ward identity 



{md + ms){md' + rris/) J dtdt MKit,t') 



E (0 I [S'lcl^iSS - DD)^cl^S' + D'^cl^{DD - SS)^cl^D'] \ 0) , (52) 

C=V,A 



which corresponds to the chiral limit of the one spin trace formalism on the lattice. Eq. 
is identical to Eq. (^) in the continuum. Since S'jclnS and D'jclnD operators cancel in 
the limit rrid = rUs, we conclude that the lefthand side must vanish as well when rrid = mg. 
The same is true for the limit md' = mg'. In the chiral limit of -EKaon —>■ 0, this behavior of 
vanishing weak matrix element could also be confirmed numerically on the lattice. 

4.2 Chiral Ward Identity of the Two Spin Trace Formahsm 

Let us define the quark propagator on the lattice as 

Gf{x,y) = {0\xf{x)Xfiy)\^) , (53) 



where / indicates the continuum flavor. In appendix A and Ref. [T^, the following lattice 
version of the continuum Ward identity Eq. (|47|) is derived 



{rrid + rris) E Gd{x, n)e{n)Gs{n, y) = e(x)G'd(x, y) + G^ix, y)e{y) , (54) 

n 

where n = {nx,ny,nz,nt) and e{n) = ^_i'jnx+ny+n^+nt _ Using Eq.(0), one obtains the 
following 

{ma + ms) (m^/ + m,/) ^ Mj[^{t, t) = Y, 

t,t' A* 

{(0 I [Xs'(7/. ® l)Xs' - Xd'(7M ® l)Xd'][Xs(7M ® l)Xs - Xd(7/. ® l)Xd] | 0) 

+ (0 I [Xs'(7m5 ® l)Xs' - Xd'(7/.5 ® l)Xd'][Xs(7M5 ® l)Xs - Xd(7M5 ® l)Xd] I 0)} (55) 

where the details of the derivation are given in appendix B and Ref. 0. Equation (^) is 
quite similar to Eq.(^). The sum Y^' in Eq. (^) indicates that in addition to the sum over 
the vector index /i we also include both one trace and two trace color contractions. From 
Eq.(p5D one can see that the right-hand side vanishes in the limit of = iris and also in 
the limit of m^i = m^i, in accordance with the continuum chiral behavior. 
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4.3 Chiral Ward Identity of the One Spin Trace Formalism 



In the one spin trace formalism, each channel requires an additional set of operators with 
different spin-flavor structure in order to realize the correct chiral behavior. From appendix 
C, we need to choose the chiral partner operator O^f^^^ partner introduced in Eq. (^Tj) as 



follows 



chiral partner 



(A X SYJ^^ + (A X S)Z:^ + (V X S)Z^^ +{VxS) 



ITR 



\1TR 
aa\bb 



with 



{A X syj^.j^^ — 
{A X sy^^f, = 

(V X 5)S - 



(V X 5)S 



1 



2Nf 
1 

1 
1 



2N 



f 



[-A ^ ^]ab;ba 



[A X 



{V X 5]^^ 



(56) 

(57) 
(58) 
(59) 
(60) 



where the notations in the righthand side of the equations are defined in Eqs. ( |18| - ^ ). Using 
the lattice Ward identity Eq. (|53), we obtain 



(m, + m,)(m,. +m,OE-^x''(^''^') =E E 



t,t 



A* {ab;ba) ,{aa;bb) 



+Xs'(7m ® I)iXsXs - XdXd)(7M ® I)Xs 



+Xd'(7M5 ® ^5){XdXd - XsXs)(7m5 ® C5)Xd' 

+Xd'(7M ® ^)(XdXd - XsXs)(7m ® ^)Xd' 

+Xs'(7m5 8) /)(XsXs - XdXd)(7M5 ® ^)Xs' 
+Xs'(7m ® 't5)(XsXs - XdXd){lt, ® ^5)Xs' 

+Xd'(7M5 ® /)(XdXd - XsXs)(7m5 ® I)Xd' 



0) 



(61) 



+Xd'(7M ® ^5)(XdXd - XsXs){'y^^ ® 6)Xd'^ 

where the details of the derivation are given in appendix C. As one can see in the above, Eq. 
(pT|) is a lattice version of Eq. (0). The key point is that in the one spin trace formalism 
a given choice of the vector and axial channel must be augmented by a particular choice of 
chiral partner channel in order to have the correct continuum chiral behavior. Therefore, 
the implementation of the continuum chiral behavior in the one spin trace formalism forces 
us to use at least eight operators rather than four operators in the two spin trace formalism. 



4.4 Chiral Limit 

In this section, we would like to discuss the leading order and the next-to-leading order chiral 
behavior of Bk- Let us consider the Fourier transform of the weak matrix elements of Bk 

C][i{E,E') = Y.^M+^Et)M^/^{t,t')exp{-iE't') (62) 
t,t' 
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C]^i{E,E') = ^exp(+^Et)Mj^(t,t')exp(-z£;'0 , (63) 



t,t' 



where M^i[^{t,t') and Mf^{t,t') are defined in Eq. (|4^) and Eq. respectively. 

For the remainder of this discussion, let us drop the label of the spin trace formahsm as 
everything apphes to both cases. Let us use the following parametrizations 0: 



Ckk^E^E') = ^^^i^^NjAKKiE,E') + Oia) (64) 

CKiE) = (ms'+md,)CKK'iE,0) = -^?^V2NfAKiE) + Oia) (65) 

Ck'{E') = {ms + md)CKK'{0,E') = -^?^V2NfAK'{E') + O{a) (66) 

E' + mj^, 

C = {ms + md)ims' +md')CKK'iO,0) = 2N]A + 0{a) . (67) 
In Ref. [|1^], it is shown that 



m 



K = \i2NffK ^ + 0{m],) , (68) 

which implies that mj^ oc + m^. Consequently, Akk' has to be expanded to 0{m'^, E"^") 

Akk'{E,E') = aoo + aioE + aoiE' + a2oE^ + auEE' + ao2E'^ + 0{E^), (69) 

where the coefficients a are functions of the quark masses up to the appropriate power. At 
this point we neglect chiral logarithms and assume analyticity of Akk'- 

For local operators in the limit of vanishing lattice spacing a — *• 0, one has 

MK{t,t')=MK{-t,-t') (70) 

and since J^Kit,t') is real it follows that Ckk' is real as well. From the definition of 
Ckk'{E, E'), we notice that 

Re[CKK'{E, E')] = function of even power of E, E' product (71) 
lm[CKK'{,E, E')] = function of odd power of E, E' product (72) 

This implies that only even powers remain which confirms that Ckk'{E, E') does not have 
a linear term in either E or E' in the continuum limit as a — [^. 

However, on a lattice of finite lattice spacing a, the operator {O2TR or ^itr general 
non-local in time. In this case, Eq. ( [fOD does not hold any more. Instead, we have the 
following 

MK{t,t) = MK{l-t,l-t) (73) 
which transfers the lattice into itself and is an exact symmetry of the non-local operator 

^Latt 
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Therefore, there exist terms with odd power in either E or E' which will be of order a. 
This implies that in Eq. (^), aio and Oqi are of order a. We finally have 

Ckk'{0,0) = ^^^Nfaoo + O{a) (74) 

= 2iVj_^ I^aoo + O{a), (75) 

which has been shown to vanish for = irid and for m^' = nidi in the previous sections. 
Therefore, the quark mass independent term in aoo must vanish and we obtain 

"00 = "00 ("^s - md){ms' - rrid') . (76) 

The 0{a) term should also vanish in the chiral limit. In the physical case of K = K' , we 
thus find for the on-shell matrix element AxKi—i'tTT'K, —if^K) 

^--hcu ^ _ ^^2^ ^ ^ 0(am;,) , (77) 

which is of the desired form. Clearly, the important condition to obtain the continuum Bk 
is that arriK <^ 1 and that the ratio ms/nid should be chosen properly (at least ms/md ^ 1) 
in the numerical simulation in order to achieve a physical value of aoo- 



5 Renormalization 



In this section, we will give a review of continuum renormalization ^\ |18|, |1^] and discuss the 
perturbative corrections to the lattice bilinear and four-fermion operators in the staggered 
fermion formulation 
Ref. 
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We adopt the powerful formulation of Ref. |]l3l. This formulation of 
was extended in Refs. 



14, PQ, 21, 22 1 to include matching all kinds of complicated 



lattice operators to the continuum operators at the one loop level. The renormalization and 

^, |2T| . The renormalization of 



matching of lattice bilinear operators was studied in Refs. 
four-fermion operators in the two spin trace formalism and their connection to the continuum 
operators was investigated in Refs. 0, |2l|, |22| . 

In addition to the above, we will also explain the renormalization of four fermion operators 
in the one spin trace formalism, including the renormalization of chiral partner operators. 



5.1 Feynman rules 



Let us briefly give the relevant Feynman rules for the staggered fermion action and lattice 
composite operators. In the general covariant gauge with a gauge parameter a (for example, 
a = 1 and corresponding to Feynman and Landau gauge respectively) , the gluon propagator 
is 



Sr.jS 



5r jA sinf^afc,,) sinf^afci. 



E/3 ^sin2(ia/c/3) 



a] 



where / and J represent the color indices in the adjoint representation of SU(3). 



(7J 
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The fermion propagator is 



. m6{p-q) + l Em sm{ap^)6{p -q + n^Ja) 
c m6{p-q) + sm{aq^)6{p-q + 7i^Ja) 



where 



TT-n.. = TC 



ir,^ — /, ^ z> , z> = unit vector along u direction (81) 

v=l 

i.e. TT^, = (0,0,0,0), 77^2 = (tt, 0,0,0), 7r^3 = (7r,7r, 0,0), TT^, = (7r,7r,7r,0), 

5{p) = (27r)^^5(p+ — n) with n = (ni,n2,n3,n4) e (82) 

n ^ 

and a and b are color indices in the fundamental representation of SU(3), a is the lattice 
spacing and — vr/a < p,q,k < n/a. 

The one-gluon and two-gluon vertices associated with the staggered fermion action are 



k 

y^{v,-<l,k)^^ = -igTl^cos{a[^+Pf,])6{p-q + k + n^Ja) 



-igT^^ cos(a[^ - q^])S{p -q + k + n^Ja) (83) 



6{p-q + ki + k2 + vr^^ /a) , (84) 



where p and g are the fermion momenta and k, ki, k2 are out-going gluon momenta. Since 
products of gauge links are present in the definition of the gauge invariant operators in 
Eg . (p!3D , the vertices of an operator will generally have external gluon lines. The operator 
vertices are 

M^^kp,-gU = ^ab4Ee^'''-^(^^®e^)Ai^e-'-^, (85) 

AB 

M^^llip,-q;k)^^ = ^agTl^^Y.^'^'^^(^i^^)ABe-^■'' 

iB,-A,)f^^^^{ak), (86) 

Ml5?i(P,-?;^i,A;2)ab = k^agrl{T^,T'},^^Y.^''''-''i'^^^)ABe--'^-^ 

- A;){B, - A,)g['^s^{ak^,ak2) , (87) 

where the superscript in curved brackets (i), i = 0, 1, 2, ■ ■ ■ denotes the number of emitted 
gluons. The function /(^^)(0 = ak) is 
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12 ,.^^.j=i 



where 



= T.<PpP^ (89) 
p 

^iin^) = Wi^, (90) 



2 

^?J(0) = l<P,f^ + ^.i^, (91) 

^?J(0) = 0-^0p/i, (92) 

C^*^) = 0-^0pA-0.^- (93) 

and = {B^ — A^)/i is defined in Eq. (|15|). All that is required to know about g'^^^-j at 
the one loop level is the following tadpole contribution 



<;b)(0>-0) = lfor/i = z/ (94) 
1 
4! 



1 i|.(A^+A.)[g + 2 J2 + 2e*^-^p^M,.^''] + h.c. for /X 7^ I. .(95) 



In order to carry out a perturbative expansion of the Landau gauge operators, we need only 
the Feynman rule in Eq. (^5)) , because all the gauge links U are replaced by the identity 
10, ^ The above Feynman rules given in analytic form are also shown graphically in 
Figure 1. [] 

Feynman rules for the four fermion operators can be easily obtained from products of two 
bilinear operators as long as the color and hypercube indices are arranged carefully. They 
will not be given here. 

5.2 Bilinear Operators 

Let us review one-loop renormalization and matching of the bilinear operators. Consider the 
bare Green's function 

{q{x)q{y)OsF) , (96) 
where Osf is a bilinear operator with spin S and with flavor F. There are 3 one loop 



Feynman diagrams [17, 18 in the continuum. We can write upto one loop 



C^fr = E (Sss'Sff^ + jfy,Z^sF-h,,)0's?^!^'^ , (97) 

S'F' \^'^) 



^ The Feynman rules given in this section are consistent with those in Ref. |T3[ |o[ |T| except for a few 
minor typing mistakes. 
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where the superscript (i), i = 0, 1, 2, ■ ■ ■ represents the number of loops and 

^SF\S'F' = ^SS'^FF'Ts log( — ) + SsS'^FF'Cg""'^ + SsS'^FF'Rs ■ (98) 

Here /i is the renormahzation scale and k is the gluon mass for the infra-red regulator. Ts 
is the anomalous dimension matrix as follows 

Ts = licTs-l) (99) 

^SlS = J^li^lulsluli, ■ (100) 

Note that the vector and axial currents have vanishing anomalous dimensions. C^""* is the 
finite constant term which depends on the regularization and renormahzation scheme. For 



the MS scheme with NDR (naive dimensional regularization) [ll7| , 

as = (4,1,0,1,4), (101) 
^cont ^ 2 10^ ^ ^^^^ ^^^^ respectively. (102) 



Rs refers to those universal terms which contain all the external momentum dependence and 
are independent of the regularization scheme. 

On the lattice, there are eight Feynman diagrams (Figure 2 (a)-(h)) contributing to the 
bilinear operator. Only half of the value of the wave-function renormahzation component 
self-energy diagrams (Figure 2 (e), (f), (g) and (h)) contribute to the lattice operator since 
the other half are absorbed in the wave-function renormahzation of the external quark fields. 
On the lattice, we have instead of Eq. ( p7|) 



2 

^SF* = "Yi^SS'^FF' + -—r:^Zsfyp)0^sfpf^ , (103) 
S'F' 

where 

Z^sfls'F' = -5ss'hF'Ts log(a«:) + C^sfIs'f' + ^ss'^ff'Rs ■ (104) 

On the lattice, is the same as in the continuum (see Eq. (PP])) since the anomalous 
dimension at the one loop level is universal, i.e. independent of the regularization and 
renormahzation schemes [|l^. Cgp!^gip, is the finite constant term unique to the lattice 
regularization scheme. For the Landau gauge definition of the axial current (75 = 7^5 and 
^F = ^5) the finite term has been calculated in Ref. [^, ^ 

^^SU.s = 12.232 . (105) 

^sf'^s'f' bilinear operators with various spin-flavor structures can be found in Ref. |20|, 
Rs is the same as that in the continuum. 

In the classical limit (i.e. at tree level as the lattice spacing a goes to zero), the lattice 
operator can be related to the continuum operator as follows 

^Latm ^ ^Conm ^ ^^^^ ^^Qg^ 
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In a perturbative calculation on the lattice, terms of order a or higher are supposed to be 
negligible as a — 0. At one loop (and higher) order the lattice and the continuum operators 
will differ and a carefully constructed mixture of lattice operators is needed to reproduce 
the desired continuum operator. From Eqs.(|97|), (|103|) , (|106|) , we can connect the lattice 
operator with the continuum operator through one loop via 



/nCont 



E 

S'F' 



SsS'SpF' + 



E 

S'F' 



(4vr) 
9 



Cont 
SF;S'F' 



yLatt 
^SF;S'F' 



/■r\Latt 
^S'F' 



Sss'^FF' + .y Sss'^FF' \n{jja) + 



(47r)2 

2 

9 ( X X /^Cont i^Latt \ 
\Oss'<JFF''^ S " ^SF;S'F' ) 



{An 



/■r\Latt 
^S'F' 



(107) 



In order to understand this operator mixing better, let us consider some important sym- 
metry properties of the bilinear operators ||2^, |1^, |14|, |22|. The distance of the lattice bilinear 
operator with spin S and flavor F ^ is defined as 



A 



,1=1 



(108) 



and corresponds to the number of links between the quark and anti-quark fields. The distance 
parity of the lattice bilinear operators is given by (—1)^. 

The massless staggered fermion action has a Ua{^) symmetry: 



X 
X 



exp[ia(75 ®^5)]X 
Xexp[m(75 (g)^5)] 



(109) 
(110) 



with 



exp[iQ;(75 0^5)] = cos(a)(/ (g) /) + ism{a){^5 
Under this transformation, bilinear operators transform as: 



Xils ® ^f)x 



cos [a] 



:-l)^sm\a)]x{l s®^F)x 
+i sin(a) cos(a) [1 + (-l)'^]x(755 ® ^F5)x- 



(111) 



(112) 



Therefore even-distance bilinear operators ^(75 ® ^f)x (i-^- A is even) can be rotated into 



XilSb ® ^F5)X , 

while odd-distance bilinear operators are invariant. This guarantees that even and odd dis- 
tance bilinear operators never mix with each other, since they belong to different irreducible 
representations respectively with respect to this Ua{^) axial rotation group. Furthermore, 
it follows that the two even-distance bilinear operators xils ® C-f)x xilsb ® ^F5)x cire 
renormalized identically (i.e. have the same anomalous dimension and the same finite term 
to all orders in the perturbative expansion). 
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We will now prove that the same is true for odd- distance bilinear operators from Ua{1] 



symmetry. This is in contrast to Sharpe and Patel, who argued in Ref. that the result 
does not follow from Ua{^) symmetry. 

Introducing a separate staggered fermion for each continuum quark flavor, we have Xu, 
Xd and Xs for three flavor QCD. Besides, for each continuum quark there are four degenerate 
flavors coming from the staggered fermion action. We then introduce a separate Ua{1) 
transformation for each continuum quark flavor 



Xu — ' exp[ia(75 ® ^5)]xu , Xu 
Xd — ^ exp[i/?(75 ® ^5)]xd , Xd 
Xs — > exp[ir7(75 O ^5)]xs , Xs ■ 



■ Xuexp[m(75 0^5)] 

■ Xdexp[i/3(75 (g)^5)] 
Xsexp[ir7(75 ® ^5)] • 



(113) 
(114) 
(115) 



Under this set of three Ua{1) transformations, the bilinear operators are transformed into 



Xs(7s ® ^F)Xd 



cos(/3) cos{r]) - (-1) sm{f3)sin{7]) Xsils ® ^F)Xd 



+i 



cos(?7) sin(/5) + (-1)^ sm{r]) cos(/3) Xs(755 ® ^F5)Xd (116) 



Even and odd distance bilinear operators still belong to different irreducible representations. 
However, now both even and odd distance bilinear operators can be rotated into 



Xs(755 ® ^F5)Xd 



This insures that for any spin-flavor structure, the two bilinear operators x(7s' ® ^f)x and 
x(755 ® ^F5)x are renormalized identically regardless of their distance. 



5.3 Four- Fermion Operators 

We now discuss the renormalization of four-fermion operators which are of the general form 



[X/i (7s ® ^f)x/2] [X/3(7s' ® ^F')X/4 



;ii7) 



where the spin-flavor structure SF; S'F' is arbitrary, allowing for 256^ = 65536 possible 
combinations. But we are interested only in those operators which contribute to the phe- 
nomenological weak matrix elements (especially Bk)- These operators are the so-called 
diagonal operators [0, ^ with S = S' and F = F' and belong to the representation / of 
the Euclidean lattice rotation group [1^, 0, |2^. For the Bk calculation, we can restrict 
the set of operators further to the following bilinear structures 



(75 ® Cf) e {(7^ ® /), (7^ ® 6), (7/.5 ® /), (7/.5 ® ^5).} 



:il8) 



'1 1^+- 

^2' 2/ 



'1 l^ 
^2' 2' 



The bilinear operators corresponding to Eq. ( |118| ) belong to J^*" 

^ respectively |2j, |13], Any linear combination of the diagonal 



and 



(I V 

V2' 2^ """" V2' 2' 

four-fermion operators with the spin-flavor structure SF; SF given in Eq. ( |118|) transforms 

as = I++ m, 11,113. 
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5.3.1 Continuum Four-Fermion Operators 

Let us consider tlie following Green's function in the continuum 

where C»g°s"* is a four-fermion operator with the spin structure S]S. From Eq. (p8|), the 
continuum four-fermion operator for (AS* = 2) is 

OlT=2 = - 75)rf"A^(l - l,)d' . (120) 



There are 10 Feynman diagrams |]T^ |18| contributing to the continuum four-fermion operator 
at the one loop level. The continuum AS* = 2 operator 0^5=2 O' 0' 01 '^^'^ t)e written 
upto one loop as 



/rxCont 



y n'Cont 



, (121) 



" (47r)2^^^=2 

where the superscript (z), z = 0, 1, 2, ■ ■ ■ represents the number of loops and 

^AS'=2 = rA5=2 log( — ) + C^T=2 + RaS=2- (122) 

K 

T^s=2 is the anomalous dimension of the operator 0^5=2 • We obtain 

rA5=2 = -4 . (123) 

^As=2 is ^ finite constant term which depends on the regularization and renormalization 
scheme. For the WS scheme with NDR fT^, 

C^S% = ^^AS=2 . (124) 

Ras=2 refers to those universal finite terms which contain all the external momentum de- 
pendence, independent of the regularization and renormalization scheme. 

5.3.2 Lattice Two Spin Trace Operators 

On the lattice the Green's function is written as 

{xixi)xiyimX2my2)OlflsF,2TR) , (125) 

where 0^'f^sf2TR belongs to the two spin trace formulation defined in section 3.1. Since 
the vector and axial channels in the two spin trace formalism satisfy the continuum chiral 
behavior separately, we need to find some linear combination of two spin trace operators 
which obey the same normalization conditions that specify the continuum operator 0^5=2- 
This is crucial so that Ras=2 in Eq- ( |122| ) cancels that of the lattice operator when we make 



a connection between the continuum and the lattice operators. 



19 



The lattice four-fermion operators for Bk in the two spin trace formahsm are 



2[xs(7,®e5)x^][x^.(7M®e5)xS'] 



N' 



^[x^(7,®e5)xS][xs'(7M®e5)x^^: 



N' 



2 [Xs (7^5 ® ^^)X'MU1,5 ® e5)xS'] 



^[x:(7/.5 ® e5)Xd][Xs'(7/.5 ® e5)Xd'] 



X pyji 



(126) 
(127) 
(128) 
(129) 



as given in Eq. (|55|-P^). They all belong to the /"'""'" representation of the Euclidean rotation 
group with discrete parity and charge conjugation symmetry [Q. There are a large number 
of four-fermion operators in the same /"'""'" representation, which will mix with the the above 
four-fermion operators. Furthermore, on the lattice color one-trace operators mix with color 
two-trace operators as in the continuum. 

Hence the general form for the one-loop renormalization of the four-fermion operators on 
the lattice is 



O 



Latt{l) 
,SF:SF 



^ij^SS' ^SS" ^FF' Op'F' 



9' 



zLatt 



(47r 



,2^ij,SF;S'F';S"F" 



^j,S'F';S"F" 



(130) 



where O^'gp.gp stands for either spin trace formulation. The subscript j of sf sf 
number of color traces and 



ryLatt 

^ij,SF;S'F';S"F" 



Ss'S"^FF'^ F' F" 



-r 



ij,SS' ln(a/t) + Rij^sS' + ^ij'!sF;S'F';S"F" 



^ij,ss' is the anomalous dimension matrix of the operator O^j'fp^gp. C^'^gp.g,p,.g„p„ represents 
the finite constant terms which depend on the regularization and renormalization scheme. 
These terms generally cause the diagonal operator to mix with operators of non-diagonal 
flavor structure at the one loop level. For B^, the total contribution of the C^-^gp.gipi.gnpi, 
mixing still belongs to the I^^ representation. R^f'Ps' represents those finite terms which 
contain all the external momentum dependence, independent of the regularization and renor- 
malization scheme as long as the anomalous dimension of the lattice operator is identical to 
that of the continuum operator as a — . 

We now have to show that there exists a particular linear combination of lattice operators 
which has the same anomalous dimension as the continuum operators at the one loop level. 

There are altogether 50 Feynman diagrams (25 diagrams each in Figure 3 and Figure 4) 
contributing to the one-loop radiative corrections to Ogpfgp2TR- The diagrams in Figure 3 
are for color two-trace operators whereas the diagrams in Figure 4 are for color one-trace 
operators. For Landau gauge operators, only 28 out of the 50 diagrams contribute to the one 
loop radiative corrections and 20 out of these 28 diagrams overlap with those of the bilinear 
operators. 

Therefore, for Landau gauge operators we only need to calculate 8 diagrams in addition 
to those diagrams for the bilinear operator: 4 diagrams from color two-trace operators are 



(131) 
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11) , (g2), (g3) and (g4) in Figure 3 and 4 diagrams from color one-trace operators are (gl), 

12) , (g3) and (g4) in Figure 4. The analytic results of the 28 diagrams for the Landau gauge 



operators are explained in appendix D and Ref. ^ . 

For Bk in the two spin trace formalism, using the scheme of Kilcup & Sharpe P, ^ |16|, 
we choose a minimum of four operators which mix with one another through the anomalous 
dimension matrix at the one loop level in order to extract the single operator that appears 
in the physical continuum 3-flavor theory 0^5=2- Defining 



O 



Bk,2TR 



{A X pyx^a 



(132) 



its anomalous dimension matrix is 



Bk 



\ {A X P)^™ ) 

/ 9 -3 -7 -3 \ 

0-62 

-7 -3 9 -3 

V -6 2 / 



(133) 



The eigenvalues of F^j^ are 
spectively 



-4, —2, 8, 16) and the corresponding eigen-operators are re- 



Ot=. 



iy X p)S 

-{v^p)Z^ + {A^p)YA 

-(^xP)^™ + (yxP)^™ 



(V^ X P)-f, + [A X PfJl + [A X P)S 



{A X PYXj^a + (^ X P)aa;ff, 
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-^iy X PflSa - iV X P)Sf. + 3(A X P)gf, + (A X P)S 



(134) 
(135) 
(136) 
(137) 



In the above, only one eigen-operator Or=-4 has the same leading logarithmic behavior as 
the continuum AS* = 2 operator upto all orders in a perturbative expansion. 

Hence by choosing (9r=-4 as the lattice version of the continuum operator, we can im- 
plement continuum leading logarithmic behavior. It is also guaranteed that the total con- 
tribution of the momentum-dependent finite terms from will cancel that of Ras=2 
when we make a connection between the lattice and the continuum operators. 

For Bx in the two spin trace formalism, the scheme-dependent finite-term matrix C^^**2r_R 
mixes Obk,2TR with a large number of other operators in the representation. In the two 
spin trace formalism, the extra operators all have different flavor structures, and so we expect 
them to be suppressed significantly in matrix elements with the external eigenstates and 
K'^ which have flavor ^5. The scheme dependent finite terms for Landau gauge operators 
Obk,2TR are explicitly 



/^Latt 
^Bk,2TR 



( 37.446 -2.9136 

28.706 

-5.25285 -2.2512 

V -4.5024 1.5008 



-5.25285 -2.2512 
-4.5024 1.5008 

37.976 -4.5043 
24.464 ) 



(138) 
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We can match the lattice operators for Bx with the continuum AS* = 2 operator by project- 
ing the eigen-operator Or=-4 out as follows 



Or°i_4 = -Pr=-4 ■ O^Br 2TR = ^^^3=2^ + terms of order a, 



^Cont(l) 



r=-4 



9' 



{An 



11. 



-iLatt 



^2 ^TB,Hfia) + -TB,-Cj^Z2TR 



o 



(1) 

Bk,2TR 



where 

Pr_4 = (l,l,l,l) 
projects out the desired eigen-operator Or=-4- 

5.3.3 Lattice One Spin Trace Operators 

Let us consider the following Green's function on the lattice 



(139) 



(140) 



{x{xi)x{yi)x{x2)x{y2)0_ 



Latt 

SF;SF,1TR. 



(141) 



where O^'^^gpu^^ belongs to the one spin trace formulation defined in section 3.2. In contrast 
to the two spin trace operators, both vector and axial channels require their chiral partner 
operators in order to respect the continuum chiral behavior on the lattice as outlined in 
section 4.3. In the one spin trace formalism, we have two questions to answer. The first issue 
is whether the chiral behavior implemented by adding a chiral-partner operator is preserved 
after renormalization. The second question is whether a linear combination of one spin trace 
operators can obey the normalization conditions which define the continuum operator. This 
guarantees that Rfj^ss' ^'i- ( |1^1| ) reproduces the corresponding continuum term when we 
make a connection between the lattice and continuum operators. The Feynman rules are 
identical to those of the two spin trace operators. Let us restrict our discussion to Landau 
gauge operators such that we only need to consider 8 diagrams (Figure 3 (gl), (g2), (g3), 
(g4) and Figure 4 (gl), (g2), (g3), (g4)) in addition to the bilinear diagrams discussed in the 
previous section (details of the calculations are given in appendix D). 

The lattice four-fermion operators for in the one spin trace formalism are given in 
Eqs. (1^), (|i^)-(^5D and their chiral partner operators in Eqs. (0), (p7D-(|60|) respectively. 
They belong to the J"*"^ representation of the Euclidean rotation group with parity and 
charge conjugation symmetry p3[. The mixing of operators is exactly the same as in the two 



spin trace formalism as described in Eq. ( 130 ). However, in the one spin trace formalism 
for Bk, there are eight instead of four operators which mix with one another through the 
anomalous dimension matrix at the one loop level 



{{V X P)llt {V X P)Zl {A X P)llt [A X P)^rj) 
and {{A X Syj^t (A x 5)^™ , (V x 5)^^, {V x 5)^) • 



(142) 



Before we go further to present the explicit results for the one-loop renormalization, let 
us consider the Ua{^) symmetry to understand the mixing of the four-fermion operators on 
the lattice better. For Landau gauge operators, these ^7a(1) rotations (Eq. ( |109|) and Eq. 



22 



( |113| , 114 , |115|) ) act on each of the bihnears of the four-fermion operators simultaneously. 
From these transformations as explained in section 5.2, we obtain the following important 



properties pi], p2[ of four-fermion operator mixing: 



(1) The distance parity of each bilinear in the four-fermion operator can not be changed 
by mixing. 

(2) The radiative corrections for the operators [xilsn ® ^F5)x][x{'1s5 ® ^F5)x] are the same 
as those for [x(7s ® ^F)x][xhs ® ^f)x]- 



(3) The renormalization coefficients of the operators [xils ® ^f)x][x{1s5 ® ^f5)x] can be 
obtained from those of [xils ® ^f)x][x{7s ^f)x] by f/A(l) rotations given in Eq. 
( [109| , |TT3|, [TT^, [TT5D. 



Numerically Ishizuka and Shizawa showed that these results are true at least at the one loop 
level for Landau gauge operators Sharpe and Patel argued in Ref. |^ that these results 
are valid to all orders in perturbation theory for Landau gauge operators. 

The property (2) in the above is useful to determine the renormalization coefficients of the 
one spin trace operators. This property insures that the anomalous dimension matrix and 
finite mixing terms of {V x P) are identical to the [A x 5") by multiplying both bilinears by 
(75 ® ^5) and the same is true for {V x S) and {AxP). Therefore, each four-fermion operator 
for Bk has the same radiative correction structure as its own chiral partner operator. Hence 
we can choose an operator basis (we call it chiral basis) such that each element satisfies 
the continuum chiral behavior and has the same radiative loop corrections simultaneously. 
In other words, the operator basis is composed of the four-fermion operator plus its chiral 
partner operator as follows 

/ {VxP + AxS)ilt \ 
{VxP + AxSyjJ, 
{AxP + VxSfXl 

V {AxP + VxS)Z!^J 

From the standpoint of group theory, this is a projection of operators to a particular repre- 
sentation of the Ua{^) group. The above basis operators {chiral basis operators) all belong 
to the identity representation of a subgroup of Ua{^) which has the following transformation 
property: 



O 



Bk,ITR 



(143) 



TT 



Xd 



exp[i-(75 ®^5)]x<i 

IT 



TT 



Xs — > exp[-z-(75 (g) ^5)]xs 



Xd 

Xs 



Xrfexp[i-(75 0^5)] 



X.exp[-i-(75 ® ^5)] 



(144) 
(145) 



The above transformation is equivalent to choosing (3 = j and 77 

(El)- 



■f in Eq. ([TT^) and 
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The anomalous dimension matrix for C_B^,iTij is the same as that for OB,i,2TR in Eq- 
( |133|) with identical eigenvalues and corresponding eigen-operators. The eigen-operator with 
eigenvalue F = —4 is 



-4,1T_R 



{A^ P + V ^ SfXl + {A^ P + V ^ S) 



aa;bb 

ITR 
aa;bb 



(146) 



Only one eigen-operator Or=-4,ir_R has the same leading logarithmic behavior as the contin- 
uum AS* = 2 operator to all orders in perturbation theory. In the continuum limit, the other 
three operators each contain at least one of the unphysical quark flavors that we introduced 
when defining the lattice theory. Therefore, these three operators can not contribute to Bk 
and can not be candidates for the continuum AS* = 2 operator. 

Hence by choosing Or=-4,iT_R as the lattice version of the continuum operator, we can im- 
plement both the continuum chiral behavior and the continuum leading logarithmic behavior 
simultaneously. It also is guaranteed that the momentum-dependent finite terms {Rfj^ss'^^ 
and Ras=2) will cancel when matching the lattice and continuum operators at the one loop 
level. 

However, the scheme-dependent finite terms C^'^\rpj^ are quite different from those in the 
two spin trace formalism. There are two classes of operators which mix with momentum- 
independent, and scale-independent, but scheme- dependent finite coefficients. The first class, 
discussed in this section, is the group of operators in Eq. ( |143| ), which are mixed by the 
anomalous dimension matrix at the one loop level. Other operators, which do not appear in 
the first class and which will turn out to be small later, make up the second class and are 
discussed in the next section. This finite mixing of the 8 operators given in Eqs. (^2]) -(^51) 
and Eqs. (|57|)-(|5(]D needs to be discussed in terms of the chiral basis operators in order to 
preserve the correct chiral behavior. Now let us present the finite mixing matrix C^'^\rpj^ in 

terms of the basis operator Obk,itr given in Eq. ( |143| ) 



/~iLatt 
^Bk,ITR 



( 36.908 -3.1442 

-0.4612 28.860 

-4.8246 -2.394 

V -4.5024 1.5008 



-4.8246 -2.394 
-4.5024 1.5008 

37.438 -4.7349 
-0.4612 24.618 / 



(147) 



where the details of calculations are explained in appendix D. Now we can connect the lattice 



operator O'^Bk itr with the continuum AS* = 2 operator (9^^™2^'' as follows 



Cont{l) 



(0) 



Pv- 
Pv- 



■ O 



(0) 

Bk,ITR 



2 

^A5=2°'' + terms of order a 



1 + 



9^ 



11. 



(47r)2 
A (P ■ Ob^^xtr 



Vbk ln(/ia) + 



^Bk^TR. 



ObkATR 



148) 



where the superscript (z), i = 0, 1, 2 ■ ■ ■ represents the number of loops, F^^ is identical to 
Eq. ( |133| ), -Pr=-4 is the projection vector given in Eq. ( |140| ) and A(P ■ Obk,itr) is the 
radiative correction which comes from the second class of operators, discussed in the next 
section. 
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5.3.4 Mixing of Operators with Different Spin-flavor Structure 



Let us discuss the mixing of operators with spin-flavor structure different from the pseudo- 
Goldstone Kaon (75 (8> ^5) in both formalisms. 

In the two spin trace formalism, there are 52 operators which mix with the original 
lattice operators in Eq. (|143|) and have flavor structure different from the pseudo-Goldstone 
Kaon. Since these operators need to be contracted with the external pseudo-Goldstone mode 
in the weak matrix elements, these weak matrix elements are proportional to a vanishing 
flavor trace in the continuum limit a ^ 0, where SU(4) flavor symmetry is supposed to be 
restored. Thus, those operators with different flavor structure {^f 7^ ^5) are suppressed by 
the vanishing flavor trace (expected to be suppressed by at least one power of a) in addition 
to the g^/i^ir)^ suppression in the two spin trace formalism. In fact the only operator that 
is not suppressed in this way is the original operator itself in Eqs. (p5|-|39|). 

In the one spin trace formalism, a group of operators mixing with the original operator 
in Eq. (|143|) which has the flavor structure different from the pseudo-Goldstone Kaon are 
suppressed only by the g'^ / {Air)"^ factor but not by the vanishing flavor trace in the contrac- 
tion with the external pseudo-Goldstone modes. This is the reason why we need to take 
into account this large group of remaining operators more carefully in the one spin trace 
formalism. 

Let us choose the chiral basis of operators described in the previous section where the basis 
operators belong to identity representation with respect to the speciflc Ua{^) axial rotation 
transformation. At the one loop level, our {{V + A) x {P + S)y^^ are mixed signiflcantly 
with the (S X + (P X {S X + (P X {V X r)i™ + {A x T)i^« 

operators which have non-negligible matrix elements with the pseudo-Goldstone Kaons: 



P ■ dB,,lTR = E ((l^ + A) X P + (1^ + A) X Sf^'^"" 



AiP-OB,„lTR 



(An) 



-1.214 (SxV + Px - 2.805 (S x V + P x 

-4.13 {SxA + Px - 2.539 {SxA + Px V)^^'^ 



-0.402 {VxT + Ax T)i°);j^ 



2.228 {VxT + Ax T)!^;^ ^ 



+ 



(149) 



where J2' means summation over two kinds of color contraction. Here we have written only 
those flnite terms corresponding to the wrong flavor operators, terms that correspond to 
A{P-dB^,iTR) inEq. O. 

The other operators (■ ■ ■ in Eq. ( |149D ) mixing at the one loop level, which are not written 
explicitly in Eq. ( |149|) , are negligible {i.e. expected to be of order a) due to the vanishing 
flavor trace in the matrix element with the pseudo-Goldstone Kaon. They are 



(T X V^)i™ + (T X , 



(150) 
(151) 
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[Xsilb X QXd Xs'{l5 X ^u)Xd' + Xs{I X ^^5)Xd Xs'{I X ^^)Xd' 
J2 [Xs(7m X ^a,y)Xd Xs'{lt^ X ^ap)Xd' + Xs{lf^ X ^au)Xd Xs'ilf. X ^<,p)Xd' 

yjXd' 

+ I] [Xs(7/.;. X ^^5)Xd Xs'{l^.u X ^p5)Xd' 
I] [Xs{li,v X ^^)xd Xs'(7a.;^ X ^p)xd' 

+ I] [Xs(7^^7x^)xd xAli^y X ^ 

1^5 jXd' 



(152) 
(153) 

(154) 

(155) 



In principle, the presence of the operators in Eq. ( |149| ) causes an extra difficulty for the 
one spin trace formalism, since there are now additional operators whose matrix elements 
must be computed. 

However, in order to evaluate the new matrix elements we can analytically express these 
one spin trace operators in Eq. ( p.49|) in terms of the small conventional set of two spin 
trace operators by Fierz transformation. Since only those operators with the same flavor 
structure (^5) as the pseudo-Goldstone mode are significant and the other operators are 
suppressed by the vanishing flavor trace in the tvjo spin trace formalism, it is possible to 
calculate approximately each one spin trace operator in Eq. ( |149| ) from the numerical results 
of the two spin trace operator measurements on the lattice. Thus the class of required matrix 
elements in the one spin trace formalism need not be expanded further 

Therefore, only those one spin trace operators whose Fierz transform has a non-vanishing 
component of either {V x p)2TR qj~ (•^ ^ py^^ are important. Now let us figure out which 
operators in Eq. (|149|) and Eqs. (|151H155"D can contribute to Bk- The Fierz-transforms of 
the operators in Eq. ( |149| ) and Eq. ( [L51| ) are 



{S X y)i^^ + {PxA) 



ITR 



-{iS-T + P)xiV + A)f^^-{{V-A)x{S-P)) 



2TR 



{S X A) 



-T + p)x{y + A)^^^ + {{y -A)x{s- p)f^^ 



{{V + A) xTf^"^ 

= -3 {{V + A)x{3S + T + 3P)f^^ 

{T X {V + A))'^'' , 

= - {{3S + T + 3P)x{V + A)f^^ . 



(156) 

(157) 

(158) 
(159) 



If we neglect those two spin trace operators with the wrong flavor structure {^p 7^ ^5)5 the 
above equations can be simplified as follows: 



(5xr)i^^ + (Px A)i™ 
{S X A)i^^ + {Px 1^)1™ 



{{V-A)xP) 
-{{V~A)xP) 



2TR 
2TR 



(160) 
(161) 
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{{V + A)xT) 
{Tx{V + A)) 



ITR 



ITR 



-^{{V + A) X P) 



2TR 







(162) 
(163) 



The Fierz transforms of the operators in Eqs. ( |151| - |155| ) contain neither {V x p)2^^ nor 



{A X p)2T"^ since those operators have non-diagonal flavor structure {i.e. the flavor matrix of 
one bilinear is different from that of the other bilinear) and they can not be Fierz-transformed 
into a desired diagonal flavor structure. Therefore, we conclude that the operators in Eq. 
( |149| ) do contain {V x p^'^'^R and [A x p)^^'^ but the operators in Eqs. ( |151| - |155D include 
neither {V x pf^^^ nor {A x Pf'^^. 

Now we can express Eq. (|149|) in terms of the non-trivial two spin trace operators as 
follows: 



A P- O 



Bk,ITR 



(47r)2 



4.122(VxP)^^«+6.418(rxP)^^« 

-1.710(A X P)2™ + 6.950(A X P)gf„ 



(164) 



The equation (|164|) can tell us how much of the radiative corrections at the one loop level 
are neglected in the leading term of the matching formula Eq. ( p.48|) . 



5.4 Tadpole Improvement 

Parisi, Lepage and Mackenzie showed that the tadpole diagrams are the main source of 
the large difference between the bare lattice coupling go{a) and the renormalized coupling 
dJTs^l^ ~ \) HH' 113' 0' ill- They suggested a mean field method for removing the dominant 
effect of tadpole diagrams. Noticing that the vacuum expectation value of the link matrix 
is smaller than 1, they proposed that the appropriate connection with the continuum gauge 
field is 



(165) 



where uq represents the mean value of the link matrix. One possible choice of a gauge- 
invariant Mo is 



Mo 



Re(-Trf/n) 



(166) 



It is the rescaled link — that should be expanded around unity. The staggered fermion 



Mo 



action can be rewritten in terms of the rescaled fields ip = a/mqX 



fermion 



'^)[x)— — '^)[x + 

Mo 



Mo 



a ijj{x)^{x) 

Mo V 



(167) 



where we take zero quark mass (M = 0) for the renormalization and the gluon action is 

1 ReTrf/n 



gluon 



MF 



(168) 
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where 



MF 



2N, 

2 ' 

9mf 



with g^F 



Re(iTrf/n) ' 



(169) 



It has been found out that ip matches the continuum quark field better than x- The pertur- 



bative expansion for Trf/n is given in Ref. p8 

1 



Re(-Trf/n) 



1 



g\ip is related to and g^-r^ perturbatively as follows 



2 

9mf 



9oiO') 



2 

9mf 



l + -9l{a)+0{gt{a)) 



A, 



l-/3o^7i,f{21n(a/i) + 21n(-^)} 



A- 



A/5 



'-,9mf + 0{giip] 



(170) 



(171) 



(172) 



X - 






X - 




y/u^X 
















Mo 



where the details of the derivation are given in Ref. |2£ 

Now let us think about the tadpole improvement for the bilinear and four-fermion opera- 
tors. For these operators, the tadpole contribution can be removed by rescaling the fermion 
fields and gauge links as follows 

(173) 
(174) 

(175) 

The tadpole improvement for the bilinear Landau gauge operators is therefore 

Off = u,OsF, (176) 
whereas for the four-fermion Landau gauge operators we have 

Ofl^S'F' = {u,fOsF,S'F' . (177) 

The only effect of tadpole improvement for the Landau gauge operators is that they are 
rescaled by a gauge invariant factor. As a consequence, the scheme-dependent finite terms 
change to 

-2.9136 

2.387 
-2.2512 
1.5008 



^Latt,MF 



:i78) 



11.207 -2.9136 -5.25285 -2.2512 

2.387 -4.5024 1.5008 

-5.25285 -2.2512 11.657 -4.5043 

V -4.5024 1.5008 -1.855 / 

/ 10.589 -3.1442 -4.8246 -2.394 

-0.4612 2.541 -4.5024 1.5008 

-4.8246 -2.394 11.119 -4.7349 

V -4.5024 1.5008 -0.4612 -1.701 / 

We would like to point out that for tadpole improvement, one has to use the renormalized 
coupling gj;jg{lJ' = ^) or gl^p rather than the lattice bare coupling g^i^a) |2^, ^ . 



C 



Latt,MF 



(179) 
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6 Summary and Conclusions 



We have presented two different metliods to transcribe tfie continuum = 2 operator onto 
tlie lattice: one is tfie two spin trace formalism (conventional method) and the other is the 
one spin trace formalism. 

The chiral behavior of the continuum AS — 2 operator needs to be respected by the 
lattice operators. Once the lattice operators are chosen, it is also important to check whether 
the lattice operators preserve the leading logarithmic behavior of the continuum = 
2 operator. Since both the lattice and the continuum operators have the same leading 
logarithmic behavior, we are guaranteed that those terms which depend on the external 
momenta cancel when we match the lattice measurement to its continuum correspondence. 

In the two spin trace formalism, it is shown that vector and axial operators satisfy the 
continuum chiral behavior separately. In the one spin trace formalism, each operator requires 
an additional chiral partner operator in order to obey the correct continuum chiral behavior. 
Hence 8 four-fermion operators are requisite for in the one spin trace formalism while 
only 4 four-fermion operators suffice for Bk in the two spin trace formalism. Furthermore, 
in the one spin trace formalism the basis operators are chosen such that they belong to the 
identity representation of the ^/a(1) symmetry group. In this representation, we can find an 
eigen-operator which has the same chiral behavior and the same leading logarithmic behavior 
as the continuum AS — 2 operator. Therefore, in both formalisms we find an eigen-operator 
which has the same leading logarithmic behavior as the continuum AS = 2 operator. 

However, the scheme- dependent finite terms differ in the two formalisms, because the 8 
operators mix with one another in the one spin trace formalism while the 4 operators in 
the two spin trace formalism. In addition, tadpole improvement through mean field theory 
is important to improve the matching between the lattice calculation and the continuum 
observable (Bk)- 

The lattice operators for Bk in the one spin trace formalism can not be Fierz transformed 
into the corresponding operators for Bk in the two spin trace formalism. However, in the 
continuum limit of a — > the contributions of the different flavor structures (for example 
((^-^^4) X 5") or {{V + A) x T^^^) vanish away. This guarantees that in the continuum 
limit a ^ 0, the matrix elements of the operators with the external pseudo-Goldstone bosons 
agree in both formalisms. 

We intend to apply both formalisms in our numerical simulations of lattice QCD. We 
will match the lattice observablcs with the continuum correspondence through the one loop 
relation given in this paper and will compare the results of both formalisms. We believe that 
this comparison will tell us how close our numerical simulation is to the continuum physics. 
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Appendix 



A Lattice Ward Identity 

Here, we derive a useful Ward identity on the lattice which will be used to explain the chiral 
behavior of four-fermion operators. This Ward identity appeared originally in Ref. [0 but 
with minor typing errors in the formula. 

The staggered fermion propagator satisfies the following relation 

[D + m]G{ni,n2) = 6n,,n, (180) 



1 



<ni\D\n2> = ^ -?7^(ni)?7(ni,n2){5„2,ni+A» - 5; 



n2,ni—fi J 



= ^Uim,n2) (181) 
G{m,n2) = <0\T[x{ni)x{n2)]\0> (182) 
Using a ^ expansion, we can obtain another representation of G{ni,n2) 
G{ni,n2) = < ni \ {D + m)~^ \ n2 > 

= -E<^ii(--n^2> 

m ^ m 

E {-ir-\-f\^\\r\Ur{n,,n2) (183) 
r(o)=ni,r(|r|)=n2 ^ ^ 

Since Uy is unitary, Uy^t2 = Uy^Uy2 is also unitary. 
Now let us consider the following Green's function 

S{ni, n2) = E Gi{ni] n)e{n)G2{n; ^2) (184) 

n 

Using Eq. (|183|) , S{ni,n2) can be simplified as follows 

S{n„n2) = EEE(-ir^""(-l)™(— )"^'+'( — )'^^'^' 
nT^T^ mi ma 

E {-ir-^{lf\Ur{n,,n2) 
r,r(o)=ni,r(|r|)=n2 

iri . . 

( _ 1 ) r(o) _ ]^ y _!_ y +1 _!_ y r| -i+i 

~^ rrii 1712 

= E(-l)"^""^(^)'"f^r(ni,n2)(-l)^(°) 
r ^ 



.{(±)|r|+i + (_i)|r|(^)|rKi^ 



nil + fn2 m2 mi 
{e(?2i)G'2(ni,n2) + G'i(ni,n2)e(n2)} (185) 



mi + m2 
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From Eq. ( |185| ), we can derive the following Ward identity 

(mi + ma) J2 Gi{ni; n)e{n)G2in; na) = e(ni)G'2(ni, ^2) + ^1(^1, ^2)6(712) (186) 

n 

Equation (|186|) will be used later for the analysis of the chiral behavior on the lattice. 

B Ward Identities for the Two Spin Trace Formalism 

We are interested in the following correlation function 

M][^{t, f) = {K{t)0,TRK'{t')) , (187) 
The relevant two spin trace operator is given by 

02TR{y) = ■^[Xs{yA){lS ® ^F)ABXd{yB)][Xs'iyc){lS ^ CF)cDXd'{yD)] ■ (188) 

Here hypercubes are labeled by capital roman letters and there is always (unless specified 
otherwise) an implicit sum over repeated indices. We define = 2y + A with A G {0, 1}'^. 
The Fourier transform Ckk' of M]^^, 

C]J^i,{E, E') = J2 eM^Et)Mf{t, t') exp{-iE't') , (189) 
t,t' 

is parameterized as in the Ref. 0. Let us evaluate 

C^fx^(0,0) = E^f^(M') (190) 

t,t' 



E [Xd(y)e(y)Xs(y)][Xs(0c)(75 ® eF)cDXd(OD)] 



y,y 



[xAOe){is ® ^F)EFX<i'{OF)][xd,{y')e{y')xAy')] 



E G^{Od; y)e{y)G,{y; Oc){ls ® ^f)cd 
y,y' 



Gd'(OF; y')e{y')GAy'\ ^e){is ® iF)EF 

1 



1 

nisi + '^d' 



G'd(0z5; Oc;)e(C) + e(D)a(OD; Oc) (75 ® ^f 



I CD 



Gd'(Oi.; OE)e(E) + e(F)G,.(0^; 0^) (75 ® ir) 



IFF 



Already at this point it is obvious that for m^ = m^ and for m^' = m^/ the above expression 
vanishes for odd distance operators, i.e. \C — D\ = \E — F\ = odd. But let us continue 



TUs + TUd 

1 

rrisi + rrid' 



GdiOo; 0c){l5ls ® Uf)cd + Gs{Od; Oc)(7575 ® ^f^^] 



CD 



Gd'iOp; Oe)(757S ® ^5^f)ef + Gs'iOF] 0e){'^S15 ® ^F^5)eF 
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1 r 



1 



, Xd' (Ob) (7575 ® ^5^f) iji^Xd' (Of) + Xs' (Of) (7575 ® -^f^s) f^Xs' (Of) 

= — ^ ^ — (a + aWcoo 

We see that the expression vanishes for = and for m^' = m^' if 

7575 = =F7575 and (191) 

This holds seperately for the {V x P) and {A x P) channel, 



(Is^^f) = (7a. ^^5) , (192) 

(75 8) ^f) = (7A.75 ® ^5) , (193) 

where the flavor is fixed to .^5 as explained in section 3.1. Notice that the vector channel 
is a distance-three operator while the axial channel is a distance-one operator. Clearly, the 
correct chiral limit is obeyed by each channel independently]^. 

C Ward Identities for the One Spin Trace Formalism 

Now we are interested in the following correlation function 

^(t, t') = {K{t)0,TRK\t')) . (194) 
The relevant one spin trace operator is given by 



OiTRiy) = ^[Xs(i/a)(7s ® ^f) ABXd'{yB)][Xs'{yc){ls ® ^F)cDXd{yD)] ■ (195) 

Now we repeat the above calculation for the one spin trace case. We have 

C]^^,{E,E') = Y,exp{tEt)Mj^'^{t,t')exp{-iE't') , (196) 
t,t' 

and 

C]Fi'{0,0) = Y.Mrit,t') (197) 



E [Xd(l/)e(2/)Xs(2/)][Xs'(0c)(75®eF)cDXd(0D)] 



y^y 



[Xs{^E){lS®^F)EFXd'{^Fmd'{yW)Xs'{y')] 



- E G'd(0z5; y)e{y)G,{y- Qe){is ® iF)cD 
y,y' 



Gd'iOp; y')e{y')GAy'; Oc){ls ® ^f)ef 

1 



1 



G'd(Oz); OF)e(^) + 6(^)^(0^; Of)|(75 ® iF)cD 
G'd'(OF; Oc)e(C) + e(F)a'(OF; Oc)}(7^®e^FF 
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Contrary to the two spin trace case, at this point it is not evident that the above expression 
vanishes for = rrid or for m^/ = nid' , because the distances between D and E and between 
F and C are not fixed. 



1 1 



TUs + rud rus' + rud' 



Ga{Od; OE)G,iOF; Oc) (757s75^ ^ 5^f^5)ef(^s ® ^f) cd + 
Gs{Od; OE)Gd'{OF; Oc)(7s ® ^f)ej.(757575 ® 6^f6)cd + 



CslOz?; OE)G's'(Oir; Oc)(7s75 ® CfC5)ef(7575 ® CfC5)cd 

1 1 



[Xd(0i?)(7575 ® e5CF)Ei.Xd'(0F)][Xd'(0c)(757.S ® e5^F)cDXd(0D)] + 



[Xd(OE ) (757575® C5^Fe5)j;jrXs'(0F)] [Xs'(0c)(75 ® CF)c z3Xd(0D)] + 

(0F)][Xd'(0c)(757s75 ibiFib) cdXs (Od)] + 



(0F)][Xs'(0c)(7575®eF6) (Od)] 

We see that the expression would vanish for = md and separately for m^/ = m^' if 

75 = (717575 = (727575 (198) 
= <J3^5^F ^ <J4^F^5 , (199) 

where cr is one of the following eight combinations 

(7= (- + + + ),(+-+ + ),(+ + - + ),(+ + +-), 

(+-—),(- + —),(--+-),( +)■ (200) 

However, this condition can not be satisfied for 75, £ {-^^ (7;ui/, ^5) ^/i5}- Instead, one 
requires an additional channel (which we call chiral partner) with different spin and flavor 
5", F' in order to achieve the correct chiral limit. The contribution is 

C'^mO,0) = (201) 



{ 



G'd(OD; 0£)G'd'(OF; Oc)(7575' 6Cf')ef(7575' CsCfOcd + 



G'd(OD; 0£)G's'(0f; 0c)(7575'75 ® ^5^F'^5) eF^S' ® ^F')cD + 



Gs{Od; OE)Gd'{OF; Oc)(75' ^ ^F')bf(7575'75 ® ^5^F'^5)cD + 
Gs{Od; Oe)Gs'{Of; Oc)(75'75 (H) Cf'C5)£;f(75'75 ® Cf'C5)cd} ■ 



We want the sum + C'j^ to vanish for = and for m^/ = rUd'- For = m^, the 
first term in the unprimed channel will have to cancel the third term in the primed channel, 
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and vice versa. Similarly, the second and fourth terms will cancel for rris' = m^/. Demanding 
that C]^j^, + C'^^ in the degenerate mass hmit, one thus obtains the following conditions 

75 = CTilbls' (202) 

75' = cr27575 (203) 

iF = ct.Uf' (204) 

^F' = cr.^F^,, (205) 

where a is one of the following combinations 

fT= (- + ++),(+-++),(+ + -+),(+ + +-), 

(+-—),(- + --),(--+-),( +)• (206) 

Let us consider o" = ( — h ++) as an example. In this case the solutions are 

(75,75') = (7m, 7/^75), (7^75,7/.) (207) 

(eF,eF') = (i,e5),(v,f^.),(e5,i) • (208) 

The other solutions can be obtained by applying the following transformations 

Ts—^Cf , Ts'-^Cf' (209) 

75 — ' -75' , ^F — > -^F' (210) 

75 ' -75' , ^F ' ^F' ■ (211) 

Especially, we are interested in the following solutions 



(75 ® Cf) = (7m ® I) and (75' ® ^F') = (7m75 ® ^s) , (212) 
(75 ® Cf) = (7m ® ^5) and (75' ®^F') = (7m75 ® /) • (213) 

such that the vector and axial channels are either distance-one or distance-three four-fermion 
operators. More generally, as in the two spin trace case, all solutions are of odd distance (i.e. 
either distance one or three). Notice that in the one spin trace formalism the correct chiral 
limit is only obeyed by the sum of both channels including their chiral partner operators, 
but not by each operator independently. 



D One-Loop Radiative Corrections 



One-loop radiative corrections for the bilinear operators are calculated in an organized way 
in several papers [|l^, |l^, ^ ^ . Here we present the important technical steps in detail 
to obtain the results for the radiative corrections for Landau gauge operators. Most of the 
diagrams for the radiative corrections for the four-fermion operators are in common with 
those for the bilinear operators. We will derive the analytic form of those diagrams which 
have nothing to do with the bilinear operators. There are eight such diagrams: (gl), (g2), 
(g3) and (g4) in both Figure 3 and Figure 4. Let us explain the technical details of calculating 
one of the diagrams (Figure 3 (gl)) in a self-contained way. Since the calculations of other 



35 



diagrams are quite similar from the standpoint of technical and mathematical difficulties. 
From the Feynman rules given in section 5.1, we can derive the naive analytic form of Figure 
4(gl) 



dp J dq J dk J dp' J dq' 
^ exp[i— ■ Aa] (75 ® ^f)ab exp[-ip ■ Ba] 

^^f A,B ^ 



la 



Eu S{p' - p + TT^n.) sin(p» 



E/3 sin^ (p'f^a) 



T4) cos(h^ - i^Ua) -6{p-^-k+^-^) 



a 



{-^9 T^y) cos[(| + (^),)a] - g + + ^) 



ta- 



il p S{q - q' + TTr^J sm{q'a) 



E/3 sin2(g/3a) 



exp[iq' ■ A'a] (75/ (g) ^f')a'B' exp[-i— ■ B'c 



N' 



f A',B' 



EA^sin^l^^; 



^sin(%^)sin(^) 
[EA^sin2(^)]2 



(214) 



where 



— _ dki rl dk2 fa dk^ dk4^ 
^ J-iL '2n J-^ ~2tx J-^ '2n J-^ '2n 



(215) 



The following technical relationship is useful to simplify Eq. 



6{p'-k + {^ + ^ + ^)) = 6ip'-k + ^ 
a a a a 



D 



with (7r^j£,'D(-l)'''' = i~'^)'''«)D'D = ® I)d' 



Da 



D 



(216) 
(217) 



This relationship is explained in detail in Ref. |Ty, |T3|, 0, Using Eq. ( |216| , ^17| ), we can 
simplify Eq. (|214|) as follows 



Gil 



I ^l >^ p f f AB ^^f A'B' 



, ki.ci 



,k^a 



cos(-^) cos(^-) sin(A;,^a) sm{kpa)F (ka) 
{6^B{ka) - (1 - a)4sin(^) sm{'^)B'{ka)) 



(-l)^'-^(7s ® ^f)ab exp{-ik ■ 5a)(-l)^"-^(7,^ ® J)^,,^ 
(7cp®/)c"C"(-1)^"'^' exp(iA; • A'a)(75' ® ^fOa'sK-I)^''^' 



(218) 
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In order to make the integral dimensionless, let us introduce a new variable such that 
= kf^a. Also we need the following identities: 



exp(ii0 ■ v4)(75 



F AB 



n^exp(i^; 



{is ® ^f)ab exp(-z0 ■ B) 



exp(-zy 



Em = 



exp(-z^) + (-l)*^^exp(+i^ 



where 



Y.EM{<P){lMS<i^^MF)AB (219) 

Y.En{-<P){isn®^fn)ab (220) 

(221) 



(222) 



It is useful to introduce the following notations for the fermion and boson propagators 

1 



Fid), m) 



Bid), m) 



J2a sin (0a) + m2 



E„4sin2(^) +m2 



(223) 



Using Eqs. (^, ^ and Eq. (|TT]) we can simphfy Eq. (PTB| ) further. 



^3 



(2vr)' 



EEEEE 

fj. u p M,N 



EM{d^)EN{-d>) cos(^^) cos(^) sin(0,) sin(0p)F2(0, 0) 



5,^5(0,0) - (1 -a)4sin(^^)sin(^|)52(0,O) 



ilSNufi ® ^Fn)cd{1CpMS' ® ^MF')c'D' 



(224) 



Let us note that this expression is divergent for M = N = i). In order to regularize the 
infra-red divergences, we need to put a vanishing mass {no) on the gluon propagator. We 
need to separate the divergent part of Eq. 



^(91) 



^3,(gl) ^3,(9l) 
^Convergent ~^ ^Divergent 



(225) 



The convergent part of G\gi) is 



^Convergent 



(27r)4 



EEEE 

fj. u p M,N 



EM{d^)EN{-d>) cos\^) sin(0,) sin(0p)F2(0, 0)5(0, ^a) 
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B{(l),0)B{(l), no) 



{iSNvix ® iFN)cD{lixpMS' ® ^MF')c'D' 



M 



[£^m(0)^m(-0) - Smo] B{(t), 0)5(0, KO) 



{iSM ® ^Fm)cd{1MS' ® ^MF')c'D' 



= 7TZn2 ^ -^ab^i'b' X! X! X! X! ^'mn{1snuij. ® ^fn)cd{1p.pMS' ® ^mf')c'D' 

\^~) I n u p M,N 

~ '^)7n2 H^ab^i'b' ^^^(^'^^ ® ^FM)cD{'yMS' iMF')c'D' ■ (226) 

y^''^) I M 

The divergent part of G^^^^ is 



jilSvn ® iF)cD{l^^vS' ® ^F')C'D' 



(1-Q;)/X^r4r^,y 1 j^,B((l),O)B((l>,Ka){^s0^F)cD(lS'®^F')c'D' 



, p, u ^ 



- {I - a) (75 (g) CF)cr>(75' ® iF')c'D' 
where is Euler number and 

i^oooo = (47r)2 /'o?xxexp(-8x) [/o(2x)]' 

JO 



(227) 



/•c 



= 4.3692- •• 



dxx 



exp(-8x)(/o(2x))' 



(47ra;)^ 



(228) 



In the continuum, we have to use the same infra-red regulator such that both infra-red 
divergences (ln(«;)^) cancel when we make a connection between the continuum and lattice 
operators. This is true as long as the anomalous dimensions at one loop are identical. The 
diagrams (g2), (g3) and (g4) in Figure 3 can be calculated analytically in the same way as 
(gl) in Figure 4. The sum of all the contributions of (gl), (g2), (g3) and (g4) in Figure 3 is 
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0000 



1) 



1 



5^ 



(47r)2 ^ 



ab-' a'b' ^mtv 

V p M,N 



■{llivNS ® ^NF - ISNvix ® ^FN)cD{lijipMS' ® ^MF' " IS'Mpp. ® ^F'm)c'D' 



■{iMS ® ^MF — ISM ® ^Fm)cd{1MS' ® ^MF' — IS'M ® ^F'm)c'D' ■ 



(229) 



Let us discuss the color one trace diagrams. The technical calculation is quite similar to 
the above procedure for the color two trace diagrams. The only thing one needs to be careful 
about is that the color structure for one color trace diagrams are different from that for two 

color trace diagrams. Performing the technical calculation as complicated as the above, the 
color one trace diagrams (gl), (g2), (g3) and (g4) in Figure 4 can be expressed analytically 
as follows 



^(gl) + ^(s2) + t^(fl3) + <^(54) 



;-ln(fi;a)^ - 7e + -^c 



0000 



p, V 



[^^'^aV'^a'b {{lixvs ® ^f)cd(75Vm ® ^f')c'D' + {isvp. ® ^f)cd(1puS' ® Cf')c'D'} 

" E^ib'-^a'b {(^t^'^S ® ^F)cD{7piyS' ® CfOc'd' + hsup O ^f)cd{1S'up ® ^F')c'D'} 



-2(1 - ln{Ka) --fE + -^0000 + 1)(75 Cf)cd{is' Cf')c'D' 

<^^<^ab"^a'b - E^ib'^a'b 



EEEE^r^'FAb"^a'b 

p, u p M,N 



(47r)2 



IpuMS ® ^Mf)cd{1S'Npp ® ^F'n)c'D' + {iSMvp ® iFM)cD{lppNS' ® CiVF')c"D'} 

ErpI rpl 
^ab'^a'b 

I 

■ {{ipvMS ® iMF)cD{'lppNS' ® ^NF')c'D' + {HSMiyp ® ^FM)cD{'yS' Npp ® ^F'n)c'D'^ 



-(i-«)7^E^. 



(47r)2 



M 



M 



[^■^'^ab'^a'b {{IMS ^ ^Mf)cd{1S'M ® ^F'm)c'D' + (75'M ® iFM)cD{lNS' ® iNF')c'D'^ 

ErpI rpl 
^ab'^a'b 
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[ilMS ® iMF)cD{.lMS' ® iMF')c'D' + {iSM ® ^FM)cDilS'M ® ^F'm)c'D']] (230) 



where 



2Nr ~ 3 



(231) 



Let us consider the diagrams in common with those of the bihnear operators. The color 
two trace diagrams (a- ■ ■), (c- ■ ■) and (d- ■ ■) in Figure 3 are identical to those of the bilinear 
operators. The color one trace diagrams (a- ■ ■), (c- • •) and (d- • ■) in Figure 4 for the Landau 
gauge operators are equivalent to those of the bilinear operators as long as the color structure 
is treated carefully. The details of the one-loop radiative correction for the bilinear operators 
are explained in Ref . [|l^, |2^, |2l| . The results are 

^2 



'^(al) 



<^(cl) 



G4 
(dl) 



9 



(4vr) 



C'i^'^ab^a'b' (^■^ ~ ~ "))(~ In(Ka) - 7ij + Fqooo + 1) 



•(7s®^f)cd(75' ®^f')c" 



D' 



+ Y^Y^Y^Y1 ^Nmill^uMSNpi, ® ^f)cd{1S' ® iF')c'D' 
fj, ly P M.N 



~(1 ~ Ci)^XM{lMSM ® C.MFm)cd{1S' ® ^F')c'D' 



M 



T^b'^a'b [(^5 - (1 - "))(- H^^^y + ^0000 + 1) 



■{is ^f)cd{is' ^Cf'! 



CD' 



+ XlfN{-f^^MSNpp ® ^f)cd(7S' ® iF')c'D' 

H u p M,N 



9 



(47r 



(1 - a)YXM{'jMSM ® iMFM)cD{lS' ® Cf')c'D' 
M 

2 r I 

a{- In(Ka)^ - 7s + -^oooo + 1) - „ ^oooo 



"C^'^ab^a'b' 



■{lS®iF)cD{lS' 


®iF') 


9^ 

(47r)'^^'^ab"^a'b 


—a{— 



'D' 



-a(-ln(Ka) - 7b + Fqooo + 1) 



^0000 



R 



R 



■{is ® ^f)cd(75' ® iF')c'D 

9' ^ ^ ^ 1 



(4vr) 



'^^'^ab^a'b' o (2 + ")^oooo(75 ® ^f)cd(75' ® ^f')c'D' 



{An 



"'^^'^ab"^a'bo (3 + a)^oooo(75 ® ^f)cd{1s' ® ^f')c'D' 



(232) 



(233) 

(234) 

(235) 
(236) 

(237) 



where 



j^B{(j),0) = (47r)2y^ rfxexp(-8x)[Jo(2x)]^ = 24.46604 (238) 
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R 



(4 



(2vr)^ 



cos(0i) (2sin^ 



2 - 2 sin2(^) + — 5(0, O)F2(0) - 5^(0, 0) 
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-5.2145- • • 



(239) 



and we follow the same notation as in Ref. 21 



The diagrams (cl) and (dl) in Figures 3 and Figure 4 were first calculated in Feynman 
gauge in Ref. |T0|. The analytic expression of the diagram (a2) can be obtained from that 
of (al) in Figure 3 and Figure 4 simply by switching the bilinear indices. The analytic 
expression of the diagrams (c2), (c3) and (c4) are identical to that of (cl) in Figure 3 and 
Figure 4. Also the analytic expression of the diagrams (d2), (d3) and (d4) are identical to 
that of (dl) in Figure 3 and Figure 4. 

Now we have the analytic expression of all the diagrams necessary for the Landau gauge 
operator. Let us introduce the following definitions 



•^(92) 



^(93) 



^(94) 



*^(9) = *^(9l) + *^{92) + <^(93) + <^(94) 



(n,c,d) 
{a,c,d) 



_ ^3 
= ^(al) 



^3 



'(a2) 



(cl) 



'(dl); 



(240) 
(241) 
(242) 
(243) 



Note that only half of the self-energy diagrams are taken into consideration since the other 
half corresponds to the the wave function renormalization of the external fields. 
Now we need the following useful SU(3) Fierz transformation identity 



1 



E^ib^ib' = -^'^ab'^a'b' + '^ab"^a'b 



(244) 



Using Eq. 



(47r)2 

-.2 



for (75 ^f) = (75' ® ^F') = ilfi ® ^5) we find 

ln(/€a)2 - 7e + Foooo + 1) [-{A X P)ab;a'b' + 3(A X P)ab';a'> 



9 



-4 



9 



(47r)2 
2 



X. 



1,22 

(0000), (0000) 



[-(AxP),,.,,,, + 3(AxP) 



ab'm'bl 



9 



(47r)2 
+4(1 - a 



(4X, 



9 



(47r 



1,11 

(nil), (nil) 
2 

^^(1111) 



+ QX, 



1,22 N 

(nil), (nil)' 



-liAx S)ab-a'b' + 1{AXS) 
D 2 



ab'\a'b 



--{A X S)ab-a'b' + T^iA X S) 

ab':a'b 

2 



+ 



(245) 



9' 



(47r)2 
-2(1 



(3 + 2«)(- ln(/€a)2 - 7^ + Foooo + 1) + SX, 
- a)X(oooo) 



1,11 I 191^^1,22 

(0000), (0000) "I" -'-^^(oooo),(oooo) 



3 1 

-(\/ X P)ab';a'fe - 2*^^ ^ P)ab;a'b' 



-12 



(4vr) 



-1,22 



|(- In(Ka) - 7e + Foooo + 1) + -^(0000), (0000) 
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7 1 

[-{A X P)ab';a'b+ X P)ab;a'> 

b 2 



9' 



(4vr) 



^^(1111), (1111) "^^(1111), (nil) 



(1 - q;)X(iiii) 



— ^(^ X S)ab';a'b — X S)ab;a'b' 



^^:^6X(nil),(llll) l^i^ ^ S)ab'-a'b - {V X S)ab-a'b'] H 



<^(a,c,d) 



where 



(47r) 



(47r)^ 



:(28.706---)(l^xP)„;.;„,;., 

: [(30.40479 • • ^(l^ x P),^,a'b - (0.56626 • • ^(T^ x P)„;.;„.6')] 



(246) 
(247) 
(248) 



X, 



1,22 

(0000), (0000) 
-1,11 



= -0.822800 • • • 
= 0.041042 • • • 



(249) 

(250) 

(251) 
(252) 



^(iiii),(iiii) 

^(nn),(iiii) = 0.023793 
X(iiii) = 0.076312 • • • 

From the above equations, we finally obtain 

(V X P)%a = 37.446(y X P)%^ - 2.9136(y x P)!^,, - 5.25285(A x P)%^ 
-2.2512(A X P)%, - 0.5381(A x S)%, - 0.2306(A x S)^, 
+0.4283(y X S)%, - 0.1428(\/ x 5)1°),, + • • • (253) 

(V X P)«,, = +28.706(\/ X P)i"j^,, - 4.5024(A x P)%^ + 1.5008(A x P)i°),, 

(254) 



0.4612(^x5)S,„ + 0.1537(^x5) 



(0) 
aa'.bb 



+ 



where the anomalous dimension terms and momentum-dependent terms are neglected, and 
the superscript represents the number of loops. The anomalous dimension matrix can be 
obtained simply by reading off the coefficients of the terms proportional to ln(«;a). 

The results for the {A x P) channel can be obtained through a procedure similar to the 
above for the {V x P) channel: 

(A X P)g,„ = -5.25285(1^ X P)%^ - 2.2512(\/ x P)^^ + 37.976(A x P)%^ 
-4.5043(A X P)^2,bb + 0-4283(A x 5)^^ " 0.1428(A x 5)i°j^,, 
-0.5381(y X S)%, - 0.2306(y x S)^^,, + ■■■ (255) 



iAxP)^I.,bb = -4.5024(yxP)l^,^,„ + 1.5008(yxP)i:^^,, + 24.464(AxP)l"^^ 



do) 



(0) 



n{0) 



bb 



-0.4612(y X 5)1°) + 0.1537(y x S)2m + 



(0) 



(256) 



where the anomalous dimension terms and momentum-dependent terms are neglected, and 
the superscript represents the number of loops. 
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The scheme-dependent finite mixing of the chiral partner operators ( {A x S)ab-ba, x 
S)aa;bb, (V X S)ab-ba and (V X S')aa;bb) JnEqs. (I^)-® Can be obtained by multiplying both 
bilinears in the above equations by (75 ® ^5). 
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E Figure Caption 



Figure 1: Feynman Rules: (a) gluon propagator, (b) fermion propagator, (c) one gluon 
vertex, (d) two gluon vertex, (e) operator with no gluon emitted, (f) operator with one 
gluon emitted, (g) operator with two gluon emitted. The dashed lines in (e), (f) and (g) 
represents gauge link between the quark and anti-quark fields as well as fiow of color indices. 

Figure 2: Feynman Diagrams of one loop correction for the bilinear operators 

Figure 3: Feynman Diagrams of one loop correction for four-fermion operators in color two 
trace form 

Figure 4: Feynman Diagrams of one loop correction for four-fermion operators in color one 
trace form 
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